home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlacon.z / dlacon
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAACCCCOOOONNNN((((3333FFFF))))                                                          DDDDLLLLAAAACCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLACON - estimate the 1-norm of a square, real matrix A
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE DLACON( N, V, X, ISGN, EST, KASE )
  13.  
  14.          INTEGER        KASE, N
  15.  
  16.          DOUBLE         PRECISION EST
  17.  
  18.          INTEGER        ISGN( * )
  19.  
  20.          DOUBLE         PRECISION V( * ), X( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DLACON estimates the 1-norm of a square, real matrix A.  Reverse
  24.      communication is used for evaluating matrix-vector products.
  25.  
  26.  
  27. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  28.      N      (input) INTEGER
  29.             The order of the matrix.  N >= 1.
  30.  
  31.      V      (workspace) DOUBLE PRECISION array, dimension (N)
  32.             On the final return, V = A*W,  where  EST = norm(V)/norm(W) (W is
  33.             not returned).
  34.  
  35.      X      (input/output) DOUBLE PRECISION array, dimension (N)
  36.             On an intermediate return, X should be overwritten by A * X,   if
  37.             KASE=1, A' * X,  if KASE=2, and DLACON must be re-called with all
  38.             the other parameters unchanged.
  39.  
  40.      ISGN   (workspace) INTEGER array, dimension (N)
  41.  
  42.      EST    (output) DOUBLE PRECISION
  43.             An estimate (a lower bound) for norm(A).
  44.  
  45.      KASE   (input/output) INTEGER
  46.             On the initial call to DLACON, KASE should be 0.  On an
  47.             intermediate return, KASE will be 1 or 2, indicating whether X
  48.             should be overwritten by A * X  or A' * X.  On the final return
  49.             from DLACON, KASE will again be 0.
  50.  
  51. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  52.      Contributed by Nick Higham, University of Manchester.
  53.      Originally named SONEST, dated March 16, 1988.
  54.  
  55.      Reference: N.J. Higham, "FORTRAN codes for estimating the one-norm of a
  56.      real or complex matrix, with applications to condition estimation", ACM
  57.      Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.